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We generalize Hardy's proof of nonlocality to the case of bipartite mixed statistical operators, 
and we exhibit a necessary condition which has to be satisfied by any given mixed state a in order 
that a local and realistic hidden variable model exists which accounts for the quantum mechanical 
predictions implied by a. Failure of this condition will imply both the impossibility of any local 
explanation of certain joint probability distributions in terms of hidden variables and the nonsep- 
, arability of the considered mixed statistical operator. Our result can be also used to determine 

' the maximum amount of noise, arising from imperfect experimental implementations of the original 

f"**) ' Hardy's proof of nonlocality, in presence of which it is still possible to put into evidence the nonlocal 

, features of certain mixed states. 
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I. INTRODUCTION 



Hardy's proof on nonlocality m has been referred to as "the best version of Bell's theorem" 0. Such a proof 
establishes, by resorting to very simple arguments which do not involve the consideration of any violation of Bell-type 
inequalities [^i 0| , a direct incompatibility between any local realistic model for almost any bipartite pure entangled 
state and the quantum mechanical predictions concerning properly chosen observables. However, the most widely 
used method to deny the existence of a local realistic model for composite states consists in the identification of a 
, Bell's inequality which is violated by the state under consideration. Contrary to the case of pure states in which any 
' entangled vector implies the violation of a precise Bell's inequality 5] , the question of which mixed states do violate a 
Bell's inequality and, as a consequence, do not admit a description in terms of a local hidden variable model is more 
complicated. In fact, within this wider scenario, there exist non-separable mixed states (that is, statistical operators 
which cannot be expressed as a convex sum of product states) which nonetheless admit a local realistic description 
and do not violate any Bell's inequality 0, ■ In brief, for mixed states the occurrence of entanglement does not in 
' general rule out a local deterministic description. Equivalently, no general procedure is known to ascertain whether 
^ . a statistical operator leads to the violation of at least a Bell's inequality H, 0, 0, 0| . 

O^' In this paper we exhibit an alternative argument not resorting to Bell's inequalities to reject the possibility of a local 
realistic description for certain mixed states. The argument is based on a reformulation and a generalization jl^J of 
"H ' Hardy's proof of nonlocality leading, via simple set theoretic arguments, to an algebraic inequality whose violation 
rS ■ by a certain mixed state implies the impossibility of a local hidden variable model for it. The novelty of the proof 
^ . derives, on one side, from the fact that it applies to a large class of mixed states, contrary to Hardy's proof which was 
restricted to pure states only, and, on the other side, that it holds for bipartite systems whose constituents belong 
to Hilbert spaces of arbitrary dimensions. The proof involves both the consideration of the trace distance of the 
considered mixed state from pure Hardy's states (that is, pure entangled states having at least two different weights 
in their Schmidt decomposition) and of the probability for precise joint measurement outcomes. The idea underlying 
our method is simply that in the vicinity (with respect to the topology induced by the trace distance) of a Hardy's 
state there exist uncountable many (non-separable) mixed states which do not admit any local realistic model. The 
usefulness of this new nonlocality argument is twofold. First, the possibility of deciding whether a given mixed state 
does exhibit genuine nonlocal features which cannot be reproduced by local classical models is extremely important 
for the theory of bipartite entangled mixed states. In fact, when this occurs one can implement efhcient quantum 
communication protocols which cannot be locally reproduced by any classical mean [l^ . Second, our result gives clear 
indications about the amount of noise which can be tolerated when performing an experimental check of nonlocality 
along the lines indicated by Hardy pi^ [l5| . More precisely, we can estimate the amount of noise (of the most general 
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kind) affecting the preparation of a pure Hardy's state so as to give a mixed statistical operator which still exhibits 
nonlocal features. 



II. THE GENERALIZED HARDY'S ARGUMENT 



Given a bipartite state vector € ® C^^, where c?i and d2 are arbitrary positive integers, let us consider its 
Schmidt decomposition in terms of appropriate orthonormal sets of states {|ai(l))} belonging to and {|/3i(2))} 
belonging to C*^ , respectively, 



IV'(1,2)) 



p^\a,{l))®\m)). 



(1) 



where the weights pi are positive real numbers satisfying the normalization condition X^i-P? = 1- Suppose that there 
exist at least two such weights which are different from each other, e.g., pi ^ P2, and let us denote any state displaying 
this property as a "Hardy state" . This is the only hypothesis which, in Hardy's proof of nonlocality , is necessary 
to exhibit a contradiction between the existence of a local hidden variable model for a Hardy's state and the quantum 
mechanical predictions for appropriate measurement outcomes. In our extension of Hardy's proof the mixed states we 
are going to consider will belong to a neighborhood of a Hardy's state, whose size will depend crucially on the values 
of pi and p2 of such a state. 

To begin with, let us consider an arbitrary bipartite statistical operator a £ B{C^'^ C^^) (that is, a positive 
semidefinite, trace class, trace-one bounded operator) whose trace distance from the Hardy state {ip) - denoted as 
D{a,\ip){'ip\) - is equal to a positive number e: 



Dia,\^){^P\) = -Tr\a-\i;){^P\\^e, 



(2) 



where \A\ = V A'^A is the positive square root of A'' A. The trace distance between two arbitrary statistical operators 
CTi and (72 represents a good measure to quantify how close are the probability distributions of any measurement 
outcome associated to the two quantum states. In fact, a well known property of the trace distance is that 



\Tr[Pai] - Tr[Pa2] \ < D{ai,a2) 



(3) 



for any projection operator P, the expression Tr[Pai] representing the probability for the occurrence of a certain 
measurement outcome when the system is associated with the state (7^. Since we already know that for any state 
vector 1^), whose Schmidt decomposition involves at least two different Schmidt coefficients, a Hardy's proof of 
nonlocality can be exhibited, our idea is that of trying to determine a neighborhood of \ip) [measured in terms of the 
trace distance D{a, \'4'){'4'\) = s] such that all mixed states a belonging to it will exhibit nonlocal features. In this 
way, we will identify a whole class of bipartite mixed states, belonging to an (arbitrary) finite dimensional Hilbert 
space, which do not admit a local realistic description and which are, as a consequence, nonseparable. 

To achieve this goal, let us recall the basic steps of Hardy's argument Q,] by resorting to the notation we used in 
our reformulation of that argument First of all, we define the following two 2x2 unitary matrices U and V 

whose entries depend on the weights pi and p2: 



U = 



1 



^/Pl + P2 



\/P2 



JP2 



V = 



1 



\/pI+pI -PiP2 



-i{ p2~P l) y/PlP2 
^/PlP2 -i{p2-pi) 



(4) 



Subsequently we consider two orthonormal bases {|a;+(l)), |x_(l))} and {|t/+(l)), |y_(l))} belonging to the two- 
dimensional linear manifold of the first subsystem spanned by the vectors {|ai(l)), |q;2(1))}, and two bases 
{|a;+(2)), |a;_(2))} and {\y+{2)), |y_(2))} for the two-dimensional linear manifold of the second subsystem spanned by 
the vectors {|/3i(2)), |/32(2))}, according to: 



l^+(l)> 
l^-(l)> 



U 



|ai(l))- 
l«2(l)) 



|y+(i)y 
|y-(i)> 



vu 



|ai(l)>' 
l«2(l)) 



(5) 



\x+{2)) 
|x_(2)) 



U 



l/3i(2)) 
1/32(2)) 



\y-m 



vu 



l/3i(2))' 
1/32(2)) 



(6) 
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The state \ip) of Eq. Q can then be expressed in three equivalent forms in terms of the basis vectors defined in 
Eqs. (|5I6|I . as: 



- zVpIpJ[|a;+(l))|x_(2)) + \x^{l))\x+{2)) ] + -Pi)|x-(l))|a;_(2)> + ^K|a,(l))|A(2)) 

i>2 

= i^pl+pl-piP2\y-{l))\x-{2))+i^^\x-imx+{2))+Y,PMmm)) 

i>2 

= i./^\x+il))\x.i2)) +i^pl+pl-p,p2\x.il))\y.{2)) +J2pMm0^i^)) ■ (7) 

i>2 

With the aim of displaying the particular set of joint probability distributions which conflict with any local hidden 
variable model, we consider the four operators Xi,Yi, X2, and Y2 having as eigenstates associated to the eigenvalues 
+1 and -1 the orthonormal vectors {|x+(l)), {|y+(l)), |y-(l)>}, {|a;+(2)), |x^(2))}, and {\y+{2)), \y^i2))}, 

respectively, while they act as the null operator in the manifolds orthogonal to the bidimensional ones corresponding 
to the nonzero eigenvalues. According to Eq. Q the quantum joint probabilities concerning the set of observables 
Xi,Yi, X2 and Y2, when the system is in the state satisfy the following relations: 

P^iXi = +1,^2 = +1) = 0, (8) 
P^{Yi = +l,X2 = -l) = 0, (9) 

F^(Xi = -i,y2 = +i) = 0, (10) 
p^(ri = +i,X2 = o) = 0, (11) 
p^(Xi = o,y2-+i) = 0, (12) 

{pi+pi-piP2r 

Since by hypothesis, the weights pi and p2 are strictly positive and different from each other [thus implying that the 
parameter a we have defined in Eq. (|13|l is strictly positive] , one is able to set up a Hardy-like proof of nonlocality T| by 
resorting to a set theoretic argument leading to a contradiction between the considered probability distributions 
of Eqs. l|S|l- H13|) and the possibility of accounting for them by means of a local realistic model where additional hidden 
variables predetermine the outcomes of any conceivable measurement. 

In order to generalize such a result to mixed states, let us consider an arbitrary statistical operator a having trace 
distance from \ip){ilj\, as defined in Eq. ((2Jl, equal to e > 0. As a consequence of Eq. (|3Jl, which gives an upper bound 
to the difference of the probability distributions associated to the different quantum states ci = ^nd 172 — cr, 

respectively, and taking into account Eqs. H8I13() . one obtains: 

P,(Xi = +l,X2 = +l)<e, (14) 

P,(yi = +l,X2 = -l)<e, (15) 

P,(Xi =-1,^2 = +1) <e, (16) 

P<,(ri =+1,^2 = 0) <e, (17) 

P<,(Xi =0,r2 = +1) <e, (18) 

P,(Yi=+l,y2 = +l)e [a-e,a + e]. (19) 

Now, suppose that there exists a more complete description of quantum systems than the one characterized by the 
simple assignment of the statistical operator a. This description is called a stochastic hidden variable model and 
it consists of (i) a set A whose elements A are called hidden variables; (ii) a normalized probability distribution 
p : A [0, 1]; (iii) a set of probability distributions Px{Ai =a,B2= &) for the measurement outcomes of any pair of 
observables Ai and B2 associated to the first and to the second subsystem respectively, defined for any value A € A, 
such that 

P,iAi=a,B2^b)= [ dXpi\)Px{Ai=a,B2 = b). (20) 

JA 

The left hand side of Eq. (|20|l gives simply the quantum probability distributions concerning the outcomes {a, b} for 
the joint measurement of the observables Ai and B2, when the system is associated with the statistical operator a. A 
deterministic hidden variable model (also known as a realistic model) is a particular instance of a stochastic model in 
which the probabilities P\ can take only the values or 1. A hidden variable model is called local if the following 
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factorizability condition holds for any conceivable joint probability distribution Px(Ai = a, B2 = b) and for any value 
of the hidden variable A S A 

Px{A, ^a,B2^b)^ Px{A, = a)Px{B2 = b), (21) 

in all cases in which the measurement processes for the observables Ai and B2 occur at spacelike separated locations. 
The locality condition imposes that no causal influence can exist between spacelike separated events. It is worth 
noticing that it has been proved 17] that deterministic and stochastic hidden variable models are completely equivalent 
when one imposes to them the locality request. For this reason, in what follows, we will deal, without any loss of 
generality, only with local realistic models reproducing the quantum mechanical predictions for the state a in terms 
of probability distributions Pa assuming only the values or 1. Finally, we will denote as the measure of any 

subset E of A with respect to the weight function p{X), i.e.. 



^(E) = / dXp{X). (22) 



To begin with, let us define the following subsets A, B, C, D of the set A of the hidden variables: 

A = {XeA\Px{Xi = 1) = 1}, (23) 

B = {XeA\Px{X2 = 1) = 1}, (24) 

C = {Ae A|PA(yi = l) = l}, (25) 

D = {Ae A|PA(r2 = l) = l}. (26) 

Suppose now that a local and realistic description exists for the mixed state a and let us consider the joint probability 
distribution Pa{Xi = 1,X2 ~ 1). With our assumptions, we have 

P^{Xi = 1,^2 - 1) = / dXp{X)Px{Xi = 1,^2 = 1) 
J A 

dXp{X)Px{X, = 1)Pa(^2 = 1) - p[A n B], (27) 



A 

where the second equality is implied by the locality condition Eq. (|21|l . and the third is a consequence of the fact 
that the product P\{Xi = 1)P\{X2 = 1) does not vanish only within the subset A O B, where it takes the value 
one, so that the whole integral gives the measure of such a set. Finally, by resorting to Eq. (|14|l . we can conclude 
that e is an upper bound for the measure of the subset AO B, that is, n P] < e. The situation becomes slightly 
more complicated when we consider Eq. H15(l and impose that there exists a local and realistic model also for such a 
probability distribution. In fact, by noticing that the only outcomes for the observable X2 are —1,0 and +1 and, as 
a consequence, that the relation P\{X2 = —1) + P\{X2 = 0) + Px{X2 = +1) = 1 holds for any A G A, we have 

P^{Yi^l,X2 = -l) - / dXpiX)PxiYi^l)Px{X2 = -l) (28) 

Ja 

- / dXp{X)Px{Yi^l)[l-Px{X2 = l)-Px{X2 = 0)] (29) 

JA 

= p[C] -p[BnC]- f dXPx{Yi = 1)Pa(^2 = 0) (30) 

J A 

Using Eqs. H15|) and (|17|l . we obtain an upper bound for the difference of the measures of the sets C and B H C 

p[C] - p[B n C] < 2e . (31) 

We can now repeat our argument for all Eqs. (|14|) - (|19|l obtaining in this way two other relations, p[D] — p[Ar\D] < 2e 
and p,[C n D] S [a — e, a + e]. Concluding, the following set of constraints on the measure of the considered subsets 
of A have to be satisfied by any local and realistic model accounting for the quantum mechanical predictions implied 
by any state a satisfying Eq. 

(32) 
(33) 
(34) 
(35) 





n[A n B] 


< 


e, 


p[C] 


-p[Br\C] 


< 


2e, 




- p\A n D\ 


< 


2e, 




p[CnD] 


e 


[a — e,a + e 
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Up to now, no constraint has been imposed on the two parameters e, quantifying the distance between an arbitrary 
mixed state a and the precise pure state |^) of Eq. and a, which specifies the non-zero probabihty of Eq. H13|l . 
Now we win show that the very assumption that a local realistic description of the implications of the state a is 
possible, implies a precise relation between such parameters. As a consequence, all states for which such a relation is 
violated do not admit any local realistic description. 

In order to find out the relation constraining the values of e and a it is useful to resort to the consideration of the 
complements A, B, C, and Z) in A of the subsets A, B, C, and D (that is, A = A — A, B = A — B,C = A — C, and 
D = A — D). We begin by taking into account Eq. and the fact that the measure of the whole set A equals 1, 
thus getting: 

^l[{A - A) U (A - B)] = ^i[A- {Ar\B)]. = l- ^[Ar\B]>l - e (36) 

In the same way, using Eq. (|35|l . we obtain: 

^l[{A - C) U (A - D)] = /i[A - (C n Z?)] = 1 - /i[C n L»] e [1 - a - e, 1 - a + e]. (37) 

Finally, by considering Eqs. 133() and 134() and resorting to simple set manipulations, we find that /i[(A — S) U (A — 
C)] - ^Ji[A -C]<2e and ^[(A - A) U (A - D)] - ^[A - D] < 2e. 

The previous relations can also be expressed in terms of the complements of the involved sets, in which case they 
take the form 

^l[A\JB] > l-e, (38) 

IJ,[B U C] - ^l[C] < 2e, (39) 

n[AuD]-fi[D] < 2e, (40) 

fi[CUD] e [1-a-e, 1-a + e]. (41) 

These new equations are more suited for deriving the desired constraint between e and a since they involve only 
the union U of subsets. In order to complete our argument, we have first of all to derive three useful relations RI-RIII. 

RI: fi[A] < n[A nD] + 2e . 

Proof. This is easily proved by noticing that 

fi[A uD]- ^i[D] = ^i[A] ~ p\A n l»] < 2e (42) 

due to Eq. ignj. ■ 

Just in the same way, taking into account Eq. (|39|l . we get 
RII: ^i\B\ < fi[B nC]+2e . 

Finally, by resorting to elementary set manipulations, we can easily prove that 

Riii: fi[A nD] + ^i[B r]C] < ^l[C yj D] + ^i[Ar] B r]C f] D]. 



Proof. It is obvious that AfMD <Z D and, similarly, that BnC C C. As a consequence {AnD)U{BnC) C (CUD) 
and fJ,[{A D D)U {B O C)] < ij,[C U D]. By the properties of any measure defined on sets, we have: 



n[{A nD)u{Bn c)] = ^l[A n d]+ fi[B nc]- iiiAn Bnc n d] < ^[c u d] (43) 

from which our conclusion holds. I 

Now we have at our disposal all the relations we need in order to derive the desired constraint between s and a. In 
fact, 

H[A UB] = ^i[A] + n[B] - n[A n < 4e + fi[A nD]+ ^i[B n C] - ^l[A n B] (44) 

< Ae + ^i[CyjD] + ^i[Ar\Br\CC]D]- ii[A n b] (45) 

< 4£ + ^[C'UL)], (46) 
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where the first majorization is implied by RI and RII, the second is implied by RIII and finally, the last inequality 
is a trivial consequence of the fact that fj.[A nBnCnD]- ^j,[A n B] < since AnBnCnD C AnB. At this point, 
using the inequalities H38|l and H41|l. we end up with the desired relation 

6£ - a > 0. (47) 

We summarize what we have just proved: 

Theorem: Let us consider an arbitrary Hardy state 2)) = J2iPi\'^ii^)) ® having two different (non-zero) 

weights pi ^ p2, and a statistical operator a such that its trace distance D{a,\ip){ip\) from the state |V')(V'I 
equals e. Then, if a local and deterministic hidden variable model exists for a, the inequality 6e — a > [where 



PlP2{Pl-P2) 

{p'i+pi-piP2y 



■] has to satisfied for any choice of the Hardy state 



The relevance of this theorem derives from the fact that it guarantees that, given a mixed state a, if there exists a 
Hardy state {ip) such that the trace distance -0(17, |V-')(V'I) is strictly less than | then no local realistic description for a 
can be given. As a consequence, any such mixed state a cannot be separable (otherwise it would admit a local realistic 
description): thus, the relation of Eq. H47|l also represents a new separability condition for any bipartite mixed state 
of the Hilbert space C^^ . 

Before concluding, let us underline two important aspects of our result. First of all, our nonlocality proof applies 
to a large class of (bipartite) mixed state. In fact, the set of pure entangled states with at least two different weights, 
i.e., the Hardy states, includes the overwhelming majority of all pure entangled states and we have just proved that 
every mixed state a belonging to an appropriate neighborhood of a Hardy state (i.e., one of size less than a/6 in the 
trace distance, where a depends on the considered Hardy state) cannot be described by a local realistic model. Thus, 
our method allows us to determine a vast class of mixed states exhibiting nonlocal features and, as a by-product, 
a vast class of non-separable states as well. Secondly, our result might be relevant from a practical point of view. 
In fact, when implementing a Hardy-experiment |14i Il5| aiming to put into evidence the truly nonlocal features of 
an entangled Hardy's state, one has to face the problem of the unavoidably imperfect preparation of such a state. 
In general one will be dealing with mixtures of states rather than with the precise pure state one wants to study. 
Actually, due to unavoidable couplings with the environment, the Hardy's state one wishes to prepare will usually 
be corrupted by different kind of noises, thus resulting in a mixed state a such as 

r7-p|^>(V| + (l-p)a, (48) 

where a is the (trace-one) statistical operator describing the noise affecting the pure state lip). Our previous argument 
shows that, notwithstanding the nonpure nature of the actual state a one is dealing with, there exists an interval of 
values for the parameter p such that one can still put into evidence a contradiction between the quantum mechanical 
predictions and those of any local realistic model. In fact, for the state of Eq. (|48|l one can repeat our earlier arguments, 
evaluate the trace distance D{(t, and determine the precise interval of values for the parameter p g (0, 1) (which 

quantifies the amount of noise in the state) for which a contradiction with locality condition is still present. 

III. CONCLUSIONS 

We have derived a generalization of the original Hardy's proof of nonlocality which works for a particular class of 
bipartite mixed states. More precisely, we have exhibited a necessary condition which has to be satisfied whenever a 
local and deterministic hidden variable model for a mixed state exists, the condition being that the state itself has 
to lay outside appropriate neighborhoods of all conceivable (entangled) Hardy's states. As a consequence, all states 
which violate the above condition do not admit any local realistic description, this in turn implying that they are not 
separable. 
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